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The magnetic vortex structure, that is present in several nanoscopic systems, is stable and can be manipulated 
through the application of a magnetic field or a spin polarized current. The size and shape of the core are 
strongly affected by the anisotropy, however, its role on the core behavior has not yet been clarified. In 
the present work we investigate the influence of a perpendicular anisotropy on the annihilation and shape of 
magnetic vortex cores in permalloy disks. We have used both micromagnetic simulations with the OOMMF 
code, and an analytical model that assumes that the shape of the core does not change during the hysteresis 
cycle, known as the rigid core model, to calculate the annihilation fields. In both cases we found that the 
annihilation fields decrease with increasing perpendicular anisotropy for almost all the structures investigated. 
The simulations show that for increasing anisotropy or dot thickness, or both, the vortex core profile changes 
its shape, becoming elongated. For every dot thickness, this change docs not depend on the dot radius, but 
on the relative distance of the core from the center of the dot. 



I. INTRODUCTION 

Among the nano- and mesoscopic magnetic structures 
that have attracted the attention of researchers in re- 
cent years stand out those that exhibit a vortex, since 
this state presents both interesting physical properties 
and a high potential for applications.^"^ Magnetic vortex 
states in nanodots are characterized by in-plane magnetic 
moments curling around a core which has magnetization 
pointing out-of-plane. Two main features are defined in 
a vortex, the circulation, i.e., the sense of the magne- 
tization curling, being —1 (+1) for clockwise (counter- 
clockwise) rotation direction, and the polarity defined 
by the direction of the core magnetization denoted by 
p = +1 (—1) for upward (downward) direction. The core 
profile nizir) (the z component of the unit magnetiza- 
tion) of a vortex in equilibrium is cylindrically symmetric, 
usually approximated by a Gaussian curve surrounded by 
a small dip (see Fig. 1).-^''^ 

A vortex configuration is the ground state of different 
nanodots with regular shape such as ellipses, squares, 
spheres, caps and disks, with lateral dimensions rang- 
ing from one hundred nanometers to a few microns, with 
some tens of nanometers thickness. ®"^^ While an external 
in-plane magnetic field that increases continuously from 
zero is applied to a disk exhibiting a magnetic vortex, 
its core will be displaced perpendicularly to the field di- 
rection, until its center reaches the disk edge. The field 
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corresponding to this limiting situation, i.e., a field that 
expels the vortex core, is known as the annihilation field. 
A further field increase will expel the vortex from the 
disk, and the saturated state will eventually be reached. 
On the other hand, when starting from a fully saturated 
state, by decreasing the field to a certain critical value 
(commonly referred in the literature as the nucleation 
field) the vortex will again be formed. The knowledge 
and control of the magnitude of these fields is a key issue 
for several applications considering the manipulation of 
magnetic vortices, such as non- volatile magnetic memory 
devices, or high- resolution magnetic field sensors. ^^"^^ 

The vortex-core nucleation and annihilation processes 
have been discussed by several authors^'^"^'^ and, in par- 
ticular, the influence of extrinsic properties on the anni- 
hilation field has been taken into account. Wu et al?^ 
investigated the role of geometrical asymmetries, find- 
ing that the annihilation of the vortex depends strongly 
on the asymmetry. The effect of the shape asymmetry 
has also been studied by Dumas et al."^^ , by measuring 
the angular dependence of the annihilation field. Miha- 
jlovic et al.^^ have shown that temperature also affects 
the reversal mechanism and the vortex annihilation field, 
while experiments by Davis et al.^'' suggest that the nu- 
cleation and annihilation fields depend on the magnetic 
field sweep rate. 

This problem was also examined from the theoretical 
point of view, within the framework proposed by Gus- 
lienko et al.^^ This model approximates the core as a 
magnetization distribution whose profile does not change 
during the reversal process. 

Some important properties of the vortices, such as the 
core size and some dynamic features, can be tailored in- 
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troducing a uniaxial perpendicular magnetic anisotropy, 
as has been recently shown^^'^^ In this the 
perpendicular anisotropy increases, important deviations 
from the vortex core profile and from the canonical mag- 
netic vortex configuration result. Beyond a critical value 
of the anisotropy (ii'^'''*), it is no longer observed a 
vortex, with the formation of a skyrmion (e.g., Fert et 
a/.^^), a structure that was found in experiments with 
BFeCoSi'^" and, more relevant to the present study, was 
also apparent in experiments with Co/Pt disks"'^'' and 
simulations. ^^'^^ 

Vortex core deformations, even under the action of a 
magnetic field, have not been so far systematically ana- 
lyzed. 

The aim of this paper is to get a better understand- 
ing of the vortex annihilation process in magnetic dots. 
For this, we have compared the description using the 
rigid vortex model with results obtained by micromag- 
netic simulations. In order to explore the effect of the 
perpendicular anisotropy on the vortex core properties, 
and concomitantly, verify the limits of validity of the rigid 
vortex model, we have introduced an anisotropic term in 
both the theory and simulations. We also characterized 
the vortex core deformations that are present in some 
simulations. 

The paper is organized as follows: after the Introduc- 
tion we describe how we perform our micromagnetic sim- 
ulations that lead us to study the annihilation fields ex- 
tracted from the hysteresis curves of disks with various 
sizes (Sec. II). Analytical calculations arc presented in 
Sec. Ill, with the inclusion of anisotropy terms into the 
rigid vortex model. The results arc contained in Sec. IV, 
and finally, in Sec. V we summarize and draw conclu- 
sions. 



II. NUMERICAL SIMULATIONS 

We investigated the hysteresis loops of individual mag- 
netic nanodots defined by their thickness L varying be- 
tween 10 and 30 nm, and diameters D from 100 to 
1000 nm, while the uniaxial perpendicular anisotropy 
Kz ranges from to 300 kJ/m^. This study was 
conducted through micromagnetic simulations^^'^^ us- 
ing the OOMMF codc.'^'^ We used a stiffness constant 
A = 13 X 10~^^J/m and a saturation magnetization 
Ms ~ 860 X 10^ A/m, the standard values used for bulk 
permalloy, taking a cell size of 5 x 5 x 5 nm^ . The maxi- 
mum anisotropics used in this work that keep the vortex 
structure are K^""" = 300, 225 and leSkJ/m^ for the 
thicknesses L = 10, 20 and 30 nm, respectively. 

For larger anisotropics, a skyrmion structure is ob- 
served, and perpendicular magnetization appears on the 
rim of the disk. For this reason, in all our calculations 
the anisotropy constant value was chosen such that the 
magnetic configuration at zero external applied field is a 
vortex configuration, as shown in Fig. 1. 

In our simulations we developed a systematic study of 
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FIG. 1. (Color Onlmc) a) profile of the vortex core corre- 



sponding to disks with D — 500 nm and L = 10 nm for 
K, = and = 300 kj/m'^; b) and c) depict the magne- 
tization for = and = 300 kj/m"^ witli L = 10 nm. 
d) profile of the vortex core with D = 500 nm and L = 30 
nm for = and — 165 kj/m^. e) and f) represent the 
magnetization for K, = and K, = 165 kj/m^ for L = 30 
nm. Note that from a) to d) the depth of the negative part 
of the magnetization (the dip) increases. 

the annihilation field that is determined from the maxi- 
mum of the derivative dM/dB in the increasing magne- 
tization branch of the hysteresis loop, which corresponds 
to the expulsion of the vortex core. All hysteresis curves 
were obtained starting from the unperturbed configura- 
tion of the disks (with the vortex core at the center), in- 
creasing the field from B = 0, in steps of AB = 0.1 mT, 
leading us to obtain the annihilation field, and finally 
reaching the magnetic saturation. In some simulations 
we observed a deformation of the vortex core. In order 
to characterize it we define 



where and ry are the sizes of the vortex core along 
the X and y axes, respectively. As shown in Fig. 1, two 
orthogonal sections (x and y directions) of the profiles 
of the vortex core passing through the core center (max- 
imum of TUz) were made. The dimensions of the core 
along the x and y directions were obtained by the full 
widths at half maximum of the repective profile fit, us- 
ing a pseudo-Voigt function. 

III. ANALYTICAL MODEL 

To obtain analytical expressions for the annihilation 
field in the magnetic nanodots we started with a model 
proposed by Guslienko et a/.^°'^^ to investigate the vortex 
behavior in submicron dots. These authors considered a 
ferromagnetic dot with a height L and a radius R that 
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(a) 




FIG. 2. Geometrical relation between the vortex core, defined 
by the dotted line, and the full dot. a) Illustration of the angle 
4>m that depends on the radius of the dot, R, radius of the 
core, h, and separation between the centers of the dot and 
core, X. b) Representation of pm, that depends on R, x, and 
the angle (p between x and pm- 



presents a vortex state with a distribution of the unit 
magnetization in cyhndrical coordinates p, (p, z given by 
TO — sin (p) + cos 6 {p) z, where^^ 



TO0 = sin 6 (p) = 



^ ( {2bp/ {b^ + p^)) p<b 



1 



p>b 



(2) 



Here 6 is the radius of the core. If we consider mag- 
netostatic, exchange, and Zeeman contributions to the 
energy, the normahzed dimensionless vortex annihilation 
field in the rigid core model proposed by Guslicnko et 
al.^^ is written as 



han R) = ^TtFi {(3) 



(3) 



where /3 = L/i?, i?o is the exchange length and Fi {(3) is 
given by 



'fit 
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A. Introducing a Perpendicular Uniaxial Anisotropy 

While the model proposed by Guslicnko et al.^^'^^ con- 
tains no anisotropy, in our calculations we include a uni- 
axial anisotropy along the z axis and focus on its effect on 
the annihilation field. We start calculating the anisotropy 
energy contribution of the system that is given by 



Wk 



LKz / (to • z) pd(j)dp , 



(5) 



where Kz > is the anisotropy constant and z is the 
easy axis. From this expression, the contribution to the 
energy due to the anisotropy comes only from the core 
region inside the dot. From Fig. 2 wc obtain 
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and 
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Using these expressions we can write Eq. 5 as 



Wk = ~2KzL 
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Wk = -KzLb^ sec 



In this expression to^ 
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and 



G represents the contributions to the anisotropy energy 
shown in the dark regions in Fig. 2a 



G = 2KzL 
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When (j)jn {x ^ R) w tt/2 or c = b/R < 1, G can 
be approximated to zero at first order of {R~ x). How- 
ever, in our calculations we considered it explicitly. If 
the anisotropy energy is normalized to M'^V, that is, 
Wk = Wk/{M^V), and using s = x/R, c = b/R and 
V = ttR^L, we obtain 



WK (s) = 
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where 17(5) = G/M'^V. We proceed by minimizing the 
magnetic anisotropy energy with respect to s and eval- 
uating in the equilibrium displacement where the vortex 
center reaches the dot perimeter. In other words, differ- 
entiating Eq. 11 with respect to s and taking the limit 
s — >■ 1, we obtain the value of the contribution of the 
anisotropy to the annihilation field 
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In this way, and adding this expression to the annihi- 
lation field given by Eq. 3, we obtain the annihilation 
field for a nanodot with perpendicular anisotropy 
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IV. RESULTS AND DISCUSSION 

Our analyses for the annihilation field are based on 
a theoretical approach and on micromagnetic simula- 
tions. From both approaches we have obtained the de- 
pendence with the perpendicular anisotropy of the anni- 
hilation fields; from the simulations the core diameters 
were also obtained. We present below the results for the 
two cases: IV A annihilation field, and IV B evolution 
of the magnetic core shape. We have observed that the 
analytical calculations result in larger annihilation fields 
as compared with the numerical simulations, however, 
there is a qualitative agreement between the results from 
both models. A qualitative agreement between theretical 
results and measured annihilation fields was previously 
report ed.^'^ 




Micromagnetic simulations 
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FIG. 3. (Color Online) Annihilation fields versus diameter 
for dots a) obtained by micromagnetic simulation, and b) 
obtained by analytical method. These graphs show the in- 
fiuence of the value of the perpendicular anisotropy on the 
annihilation field for difllerent diameters. 



A. Annihilation fields 

The values of the annihilation fields Ban as a function 
of disk diameter, for different anisotropies Kz, obtained 
from both methods are given in Figs. 3, 4 and 5 for 
L = 10, 20 and 30 nm, respectively. 

The two sets of results agree in the fact that, as the 
diameters of the disks increase, the values become less 
dependent on the anisotropy (Figs. 3, 4 and 5). 
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FIG. 4. (Color Online) Annihilation fields versus diameter for 
dots a) obtained by numerical calculation and b) obtained by 
the analytical method. 

Increasing the anisotropy, which results in larger core 
sizes, one is led to lower annihilation fields. 

For small diameters, the effect of the anisotropy is more 
noticeable than for the large disks, both in the theory 
and simulations. In the simulations, the sensitivity of 
the annihilation field to changes in anisotropy increases 
with the value of Kz- 

The disagreement between the two methods can be 
related to a deformation of the core observed in the mi- 
cromagnetic simulations. 
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FIG. 5. (Color Online) Annihilation fields versus diameter for 
dots with a) and obtained by numerical calculation and b) 
obtained by the analytical method. 



B. Evolution of the magnetic core shape 

In the search of the effect that gives rise to the dis- 
agreement between the values of Ban derived by ana- 
lytical method and by numerical calculation, as we men- 
tioned above, we considered the possibility that this effect 
would also cause a deformation of the vortex cores. Our 
micromagnetic simulations show that this deformation is 
in fact present, as shown in Figs. 6 and 7. Of course, 
the rigid vortex model breaks down in the cases where 
important deformations of the vortex core are observed. 

We have therefore analyzed the evolution of the shape 
of the magnetic vortex core as it moves towards the edge 
of the disks, under the influence of an applied magnetic 
field; the effect of varying the value of a perpendicular 
anisotropy Kz was also investigated. The simulation re- 
sults show a gradual deformation of the vortex cores, that 
can be quantified; the vortex cores change from a cir- 
cular shape to a nearly elliptical ("banana-like") shape. 
Figs. 6 and 7 show images of the vortex core and the 
core profile for L = 10 nm and L = 30 nm; the change in 



shape of the vortex core is very clear. 

Comparing Fig. 1 and Figs. 6 and 7, we observe that 
the vortex core deformation increases as the core moves 
away from the center. The deformation is zero (i.e., the 
core is circular) for B = and maximum when the vortex 
core is on the edge of the disk, in the instant immediately 
before the vortex annihilation. 

To clarify this point we have drawn the magnetization 
profile of the core for fields close to the annihilation field 
for dots of L = 10 and 30 nm, as shown in Figs. 6a 
and 7a, respectively. The dark region (blue online) in 
these figures represents the core region. The profile of the 
vortex core at the center of the disk is shown in Fig. 1. 
The Figs, la, lb and Ic show the effect of the anisotropy 
on a dot of diameter 500 nm with L = 10 nm, and Figs. 
Id, le and If for diameter 500 nm, with L = 30 nm. 

For L = 10 nm the core keeps a nearly circular 
shape, the core deformation reaches about 10 % for zero 
anisotropy {Kz = 0), and for Kz = 300kJ/m'^ the defor- 
mation is around 30%. Comparing Figs, la, lb and Ic 
with Fig. 6, it is evident the vortex core deformation in 
the latter, due to an increase of its size along the y axis 
and the depth of the magnetization dip. 

However, for L = 30 nm, the smallest core deformation 
is 30%, that corresponds to anisotropy Kz = 0, and a 
larger core deformation is observed for Kz — 165kJ/m"^, 
of about 100%. Therefore, the core loses its circular 
shape for large anisotropy constant values Kz, leading 
to a core of roughly elliptical section that is not well de- 
scribed by the rigid vortex model. This is shown in Figs. 
6 and 7, where it is more evident the deformation of the 
vortex core, as well as the variation in the magnetiza- 
tion dip. It is important to note that a higher anisotropy 
constant for L = 30 nm will result in a magnetic configu- 
ration that is not a vortex, as observed above. Therefore 
Kz = 165kJ/m'^ is the largest value that can be consid- 
ered in order to obtain a vortex state. 

We observe that the core deformation is not present for 
B = 0, as shown in Fig. 1. However, as the field begins 
to increase, the process of deformation of the core sets 
in, as can be seen in Figs. 8 and 9, plotted from values 
obtained using Eq. (1); they are shown as a function of 
relative core positions {pcore/ R), defined as the ratio of 
the distance of the core center from the center of the disk 
(Pcore) divided by the disk radius (i?). 

Figures 8 and 9 show different stages of deformation. 
In the case where the vortex core is located below 25% of 
the radius of the disk, the deformation can be neglected. 
From this point onwards the deformation begins to in- 
crease, and the maximum is reached when the core ap- 
proaches the edge of the disk. Note in Figs. 8 and 9 that 
for each thickness (10, 20 and 30 nm) four curves were 
plotted for different diameters (400, 500, 750 and 1000 
nm) and that these curves overlap; it appears to exist a 
scaling law for the deformation of the core, i.e., for each 
disk thickness, the deformation does not depend on the 
disk diameter, it depends only on the relative position of 
the vortex core. 
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FIG. 6. (Color Online) Change of the shape of the vortex core 
for L = 10 nm, D — 500 nm and K~ = 300 kj/m'^ immediately 
before the annihilation, a) profile of the core along the x axis 
(red dotted line) and along the y axis (blue continuous line), 
b) representation of the disk, c) a detail of the disk section 
close to the vortex core. 
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FIG. 7. (Color online) Change of the shape of the vortex 
core for L = 30 nm, D — 500 nm and Kz = 165kJ/m^ at 
positions immediately before annihilation. In a) profile of 
the core along the x axis (red dotted line) and along the y 
axis (blue continuous line), in b) image of the disk and c) a 
detail of the disk. 



V. CONCLUSIONS 

In summary, by means of an analytical model and nu- 
merical simulations we have obtained the annihiliation 
fields for dots of different heights and anisotropy con- 
stants. In all cases, the annihilation fields decrease with 
increasing anisotropy constant and with increasing 
disk diameter. The values of and disk height L have 
a stronger effect on the annihilation fields of the smaller 
disks. However, the anisotropy and the thickness operate 
in an inverse way on the annihilation field; whereas the 
increase in height increases the annihilation field, from 
the analytical results using the rigid vortex model, the 
increase in anisotropy decreases this field. Finally, we 
have shown the variation in the deformation of the vor- 
tex core (5 as a function of the perpendicular anisotropy 
and disk height. The occurrence of this deformation ev- 
idenced in the micromagnetie simulations suggests that 
it has to be taken into account in the description of the 
dynamics of the magnetic vortices. The deformation of 



the core does not scale with the radius of the disks, it is 
only related to the relative position of the core. 
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